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We consider the dynamics of atomic and field coherent states in the non-resonant Dicke model. At 
weak coupling an initial product state evolves into a superposition of multiple field coherent states 
that are correlated with the atomic configuration. This process is accompanied by the buildup and 
decay of atom- field entanglement and leads to the periodic collapse and revival of Rabi oscillations. 
We provide a perturbative derivation of the underlying dynamical mechanism that complements the 
rotating wave approximation at resonance. The identification of two different time scales explains 
how the dynamical signatures depend on the sign of detuning between the atomic and field frequency, 
and predicts the generation of either atomic or field cat states in the two opposite cases. We finally 
discuss the restrictions that the buildup of atom-field entanglement during the collapse of Rabi 
oscillations imposes on the validity of semi-classical approximations that neglect entanglement. 
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I. INTRODUCTION 

In cavity quantum electrodynamics pQ , where the con- 
finement of atoms results in coherent coupling to a sin- 
gle field mode, optical signatures such as non-Poissonian 
photon statistics [2] or vacuum Rabi oscillations [3 allow 
for the direct observation of quantum effects on light- 
matter interaction. A fundamental consequence of field 
quantization is the collapse and revival (CR) of Rabi 
oscillations in a resonant cavity [4]. This effect, de- 
scribed by the Jaynes-Cummings model [5 , involves the 
generation of atom-field entanglement and non-classical 
"Schrodinger cat" states of the photon field [6]. In non- 
resonant cavities, on the other hand, the preparation of 
field cat states relies on linear frequency shifts induced by 
the atom- field coupling [7 . Additional effects arise in sit- 
uations beyond weak coupling or resonance [8], but also 
for the superradiant quantum phase transition in atomic 
ensembles [9] which has a close connection to quantum 
chaos and critical entanglement [TO] . 

In this paper we analyze the CR dynamics of atomic 
ensembles in the non-resonant Dicke model [TT] . We iden- 
tify the relevant CR mechanism that results from the 
weak-coupling correction to the bare atomic and field 
frequency. In difference to the behavior of a single atom 
studied in the Jaynes-Cummings model, more complex 
CR patterns are expected for atomic ensembles [12]. So 
far, they have been discussed only in the rotating wave 
approximation (RWA) p~3j[14], which is restricted to the 
near-resonant case. 

The analysis of the non-resonant case provided here 
shows the importance of two different time scales for the 
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CR dynamics. They are associated with the dynamical 
splitting of either field or atomic coherent states, and 
the subsequent generation of field or atomic Schrodinger 
cat states. The effectiveness of the different dynamical 
mechanisms depends on the sign of detuning: Atomic 
(field) cat states are generated predominantly if the field 
frequency is larger (smaller) than the atomic transition 
frequency. The former (latter) situation is further char- 
acterized by the buildup of significant atomic squeez- 
ing (atom-field entanglement) during the initial collapse 
phase. In both cases a periodic CR pattern develops on 
long time scales. 

To understand these effects we proceed as follows. We 
first describe, in Sec. [TTJ the principal behavior using nu- 
merical results for atomic expectation values, the entan- 
glement entropy and phase space distribution functions. 
These results establish Rabi oscillations and CR dynam- 
ics for the non-resonant case and indicate the evolution 
of an initial product state into a quantum superposition 
with large atom-field entanglement. We then deduce this 
behavior from the non-resonant weak-coupling pertur- 
bation theory developed in Sec. Ill and App. [A] as the 
equivalent to the RWA at resonance. Perturbation the- 
ory allows for a clear identification of the relevant mech- 
anisms and predicts the characteristic structure of the 
wave function as a quantum superposition of multiple 
classical field states. In Sec. [IV] we return to numerical 
calculations for the opposite cases of small or large field 
frequency. According to the two different time scales that 
we found in perturbation theory we will observe the re- 
alization of "Schrodinger cat" -like states of the atom or 
the field. Since atom-field entanglement plays a decisive 
role in the CR dynamics, we discuss in Sec. [V] the conse- 
quences for the standard semi-classical approximation to 
the Dicke model, which is found to encounter problems 
even close to the classical field limit. Our findings are 
summarized in Sec. I VII 
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FIG. 1. (Color online) CR dynamics of three atoms (j — 3/2) 
for a/A = 0.01, A/A = 0.01, with O = tt/2, a = 5.0 in 
the initial state at t = 0. Upper row and lower left panel: 
Atomic expectation values (J x (t)) (black curve) , (J z (t)) (blue 
curve) over different time scales, covering three orders of mag- 
nitude. Panels (b)-(c) show envelope functions because the 
fast atomic oscillations visible in panel (a) can not be resolved 
in the pictures. Lower right panel: Entanglement entropy S(t) 
(red curve) and field variance A/(t) (green curve). 



II. COLLAPSE AND REVIVAL DYNAMICS 

Atomic ensembles are described in the Dicke model [11] 

H = -AJ Z + A(a f + a)J x + £la ] a , (1) 

as a pseudo-spin of length j (using angular momentum 
operators J X1 J z \ which represents an ensemble of 2j 
two-level atoms with transition frequency A coupled to 
a bosonic field mode with frequency Q (using ladder op- 
erators a^)). We consider this model in the two different 
non-resonant cases 0<A and £1 ^> A. 

An example of CR dynamics for Q <C A is given in 
Fig. [T] In this and all following examples, the system at 
time t = is prepared in the product state 



Wo)> = |0 O > ® N> 

of an atomic coherent state 



(2) 
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and a field coherent state 

|ao) = exp[ao^ — <^o a ]|0) 

oo n 
n=0 



(4) 



using the J^-eigenstates | j, m) and the a^a-eigenstates \n) 
(cf. Ref. [15 J. In the initial state it is (J z ) = jcos#o? 
(Jx) — jsin^o and (a) = ctQ. We assume that |ao| ^> 1 
for a classical field, and generally choose G M, > 0. 

Fig. [I] shows the atomic observables (J x / Z (t)) over dif- 
ferent time scales. They have been calculated from nu- 
merical time-propagation of the wave function using the 



Chebyshev technique [16 j. Up to 10 3 bosons have been 
kept in the calculations to prevent errors from the trunca- 
tion of the bosonic part of the Hilbert space. All numer- 
ical data shown are exact in the sense that the relative 
error is on the level of machine precision. 

On a short time scale (panel (a)), we observe fast os- 
cillations with the atomic frequency A. The amplitude 
of the (J^t)) -oscillations is of order A. On a longer time 
scale (panel (b)), the atomic oscillations in (J z (t)) are 
modulated by oscillations with frequency ft (note that 
we show envelope functions whenever the fast atomic os- 
cillations are not resolved in the pictures). We call them 
Rabi oscillations in analogy to the resonant case because 
they arise from the coupling of the atoms to a classical 
field. In the non-resonant case 1]«A they appear with 
the field frequency Q. The collapse of Rabi oscillations is 
observable over the first 5-10 field periods 2it/Vt (panel 
(c)), before they reappear on an even longer time scale 
(panel (d)), with a revival time of T r /{2k/&) w 100 field 
oscillations. In contrast to the resonant case, a periodic 
CR pattern of Rabi oscillations evolves. 

The CRs are accompanied by the periodic buildup and 
decay of atom-field entanglement (lower right panel in 
Fig. [I]) , which we measure through the entanglement en- 
tropy 



-Tr[p r lnp r ] . 



(5) 



It is obtained from either the reduced atomic or field den- 
sity matrix p r , which both give the same value according 
to the Schmidt decomposition [17] . From the initial prod- 
uct state with 5 = 0, entanglement is generated in the 
collapse phase until S is close to the maximal possible 
value ln(2j + 1). The revivals coincide with entangle- 
ment decay, as the wave function returns to a product 
form. A similar behavior is found for the field variance, 
which we define as the product 



A / = i(A[at + a]A[i(at-a)] 



1/2 



(6) 



of the uncertainties of the field operators (a"*" + a) / y/2, 
i(a t -a)/V / 2, where AA = (A 2 ) -(A) 2 as usual. In Fig.[l] 
it signals that the initial coherent state with minimal 
Af = 1/2 evolves into a field state with large variance in 
the collapse phase. 

To identify the structure of the bosonic field we show 
in Fig. [2] the field Husimi function [15] 



Q f (a) 



(7) 



By definition, it is the probability of finding the field in 
the coherent state \a). For A = 0, points in the phase 
space rotate on circles around the origin a = 0. We ob- 
serve that the Husimi function is a superposition of sev- 
eral Gaussian peaks indicating the splitting of the initial 
into multiple coherent field states in the collapse phase. 
Partial revivals, e.g. at t = Tr/2 ~ 50 x (27r)/fi, occur 
when some but not all of the coherent field states merge. 
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FIG. 2. (Color online) Field Husimi function Q/(a) in the 
collapse phase for three atoms with parameters from Fig. ^ 
at times t/ (27r/0) = 10, 25, 50 corresponding to the arrows in 
the lower left panel in Fig. [I] The dashed circles have radius 
ao. The color encoding runs from white for Q = through 
blue to red for the maximal Q- value in the respective picture, 
as shown in the color bar on the right. 



At a full revival t = Tr (not shown), the Husimi function 
consists again of a single Gaussian peak at which in- 
dicates the revival of the initial state. This explains the 
behavior of the field variance Af. 

We observe here the first of two CR mechanisms. An 
initial coherent field state splits because the rotation fre- 
quency in oscillator phase space depends on the atomic 
configuration. The field-induced collapse phase coincides 
with maximal entanglement between the atomic and sev- 
eral coherent field states whose phase space position dif- 
fers by a finite angle. The large entanglement indicates 
an incoherent superposition rather than a field cat state 
that requires a coherent superposition [18 j. Periodic re- 
vivals occur when the relative angles approach zero again. 



III. NON-RESONANT PERTURBATION 
THEORY 

After the discussion of the first example, we now derive 
the CR mechanisms using perturbation theory. The cen- 
tral result will be that, under the appropriate conditions 
for the perturbative treatment specified below, the wave 
function has the structure 



(8) 



where \a m (t)) are a set of atomic (i.e. spin) states that 
occur together with field coherent states \a m (t)}. We 
note that the field remains essentially classical during 
time evolution since only coherent states occur, but it be- 
comes entangled with the atomic ensemble because of the 
dependence of a m (t) on m: Each classical field state is 
"tagged" by the associated atomic configuration \cr m (t)). 



A. Derivation of the perturbative wave function 

The unperturbed eigenstates at zero coupling (A = 0) 
are the product states \m) <g> \n) of J^-eigenstates \m) 



and field Fock states |n), with energy E^ nn = -mA + nl]. 
From standard non-degenerate second order perturbation 
theory we obtain the energy correction as 



E {2) 



A' 



■E 

jJLV=±l 



(j — /jLm)(j + \im + l)(2n + v + 1) 
fiA — 



(9) 

where the sum contains contributions from the four states 
|ra±l) (g) |n± 1) contributing in second order through the 
interaction term \(aJ-\-a) J x . The correction to the eigen- 
states is given by a similar expression, but we need to 
keep only the leading first order terms in A (further per- 
turbative results are given in App. |A|). Both expressions 
can be combined into a unitary time-evolution operator 
that gives the perturbative wave function as 



\m) 



(2) 



17+ 



exp[-Lfft]tf|^(0)) , 



with an effective Hamiltonian 

H = -AJ Z + Vta ] a - u E (2a t a + l)J z - u s J z 
and a unitary transformation of states 



(10) 



(11) 



U = exp 



(o T - a)J x 



A v '— A 
We here introduced the two frequencies 

A 2 A A 2 fi 



i^(<f+a)J v . (12) 



0J E 



2 (a 2 - n 2 ) ' 



UJ S 



2(A 2 - n 2 ) 



(13) 



that appear in H as a consequence of the atom-field cou- 
pling, and dropped a constant term + 1). Due 
to the unitary form of the perturbative result it re- 
mains valid for long times t and large |ao|, provided that 
X\a \ < |A 2 -^ 2 |. 

The central information about the non-resonant CR 
mechanism is contained in the two time scales 



\oj e \ 



7T 
1^1 



(14) 



in the effective Hamiltonian H. The "entangling" time 
Te is associated with the term a^aJ Z: which gives an 
energy correction ex ran and is the origin of the atom- 
dependent field splitting observed in Fig. [2} The "squeez- 
ing" time Ts occurs with J 2 . Since this term affects only 
the atomic ensemble, no additional entanglement is gen- 
erated. Instead, it leads to the squeezing of atomic co- 
herent states and splitting into atomic cat states [13j [19] . 
The ratio Te/Ts = Q/A determines which term domi- 
nates the initial dynamics. A similar perturbation theory 
for a simplified model in RWA does not distinguish be- 
tween the different time scales [14 (see also App. [A]). 

We can now construct the perturbative wave function 
in the form of Eq. (J8|, starting from the initial state 
1^(0)) = \0o) \ao) used throughout this paper. Under 
the assumption \cto\ ^> 1 we can replace the operators 

+ a, i(aJ — a) in the unitary transformation U from 
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Eq. (12) by the scalars 2Reao, 2Imao respectively. The 



error of this replacement is of order l/|ao|. Then, U 
reduces to a spin rotation operator of the form 

i?(a, b) = exp i(a J x - bJ y ) (a, b G R) , (15) 

and we get 

'4cJs-Ima 4cj£;Rea^ 



|cr)(g)|a) 



R 



A 



A 



)k>]®|a) (16) 



for every atomic state \a) and a field coherent state with 
\a\ > 1. 

We note that general atomic coherent states can be 
defined through 



which gives 



R(9 sin 0, 9 cos j) 



(0,0|J*|0,0)=jcos0, 

(0, (j)\J x \9, (j)) = j sin # cos . 



(17) 



(18) 



Atomic coherent states remain coherent states under ro- 
tation. In particular for the initial state we have 



U 



\0 ) ® \ao) = \e + 5e) ® |a ) 



(19) 



from the relation R(0, 0')\6) = \0 + 0'), a simple rotation 
of the atomic coherent state by the angle 



69 = 



2\Aa 

a 2 - ft 2 



(20) 



Here, we still assume for simplicity that ao G M. 

Since the effective Hamiltonian H is diagonal in the 
J^-eigenstates |ra), application of exp[— L^Tt] to the state 
in Eq. (19) rotates the field component of the different 
J^-contributions. The operator a^a generates a rota- 
tion of field coherent states of the form exp[i^a]|a) = 
|exp(i£)a), such that we have 

exp[— \Ht}\m) (g) \a ) 

= e i<-(A+^)+m 2 o; s ) ^ ^ ^(ft) , (21) 

with 

a m (t)=a e- it ^- 2rnuJE) . (22) 
We note the scalar product 

\(a m (t)\a m >(t))\ = exp |a | 2 (l - cosfo mm /)) (23) 



between two field coherent states enclosing the finite 
phase space angle 5a mrn > = 2i\(m — m / )(t/Tg;). 

For each of the states in Eq. (21) the inverse transfor- 
mation leads again to a spin rotation as in Eq. (16), 
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FIG. 3. (Color online) Perturbative result for the CR dynam- 
ics of three atoms with parameters from Fig. [I] Shown are 
the atomic expectation values (J x (t)), (Jz(t)) obtained with 
Eqs. (j8|-(25), corresponding to panels (b)-(d) in Fig. [I] 



but the arguments of R(-,-) now depend on a m (t). We 
have W\m) <8> \a m (t)) = \cr m (t)) <g> \a m (t)) with 

/ 4cJ5 4c<j£; \ 

Wm(t)} = R[ — Ima m (t),-— —Rea m (t)j\m) . 

(24) 

Collecting all results, we finally see that the pertur- 
bative wave function |V^))(2) nas indeed the structure 
proposed in Eq. (|8|, with individual terms given by 
Eqs. p2|, ([24]) and 



The coefficients (m\9) of an atomic coherent state are 
given in Eq. (J3j) . We note that the spin states \cr m (t)) are 
not orthogonal since the rotation with •) depends on 
a m (t). This effect is of order A and not an artifact of 
perturbation theory. 



B. Collapse and revivals in perturbation theory 

Explicit expressions for the values of atomic and field 
observables can be derived from the above equations, but 
we do not show them here since they are rather lengthy 
and not very illuminating. Instead, let us focus on the CR 
dynamics in the situation 0«A addressed in Figs.[T][2j 
The perturbative result for atomic expectation values is 
shown in Fig. [3] Comparison with Fig. [I] shows that 
perturbation theory is in excellent agreement with the 
numerical data. 

In the initial dynamics, for t <C T^T^, differences 
between the coherent state parameters a m (t) are neg- 
ligible. Therefore, the wave function has product form 
\ip(t)) « \a(t)) <g> |a e" im ). The atomic state \a(t)) 
is obtained from the initial atomic coherent state |0q) 
through three rotations around different axes: The first 
rotation from Eq. (19), the second rotation from the ef- 
fective Hamiltonian as exp [it A J^], and the third rotation 
from Eq. p4| ) which presently does not depend on m. 
Therefore, \o\t)) is itself an atomic coherent state as in 
Eq. (17). From Eq. (18) we obtain the atomic expecta- 



tion values as 



(J z (t)) = j(cosAtsin((S0cosfit)sin(0 o + 59) 

+ cos(S9 cos ttt) cos(<9 + £0)) , (26) 
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(Jx(t)) = j(cos At cos(S9 cos fit) sm(9 + 50) 

- sin(50 cos fit) cos((9 + 50)) . (27) 



These expressions describe, through the term 50 cos fit, 
Rabi oscillations with frequency ft. Their origin within 
perturbation theory is the dependence of the final rota- 
tion in Eq. ([24j) on a m (t) « aoe~ iQt for t<^T E . 

Since 1]« A, the term a^aJ z determines the CR dy- 
namics for times t ~ Tg C T^. The splitting of the 
initial coherent field state \ao) into the 2j + 1 coherent 
states \a m ) is the source of entanglement with the atomic 
J^-eigenstates \m). Because different \m) -states are or- 
thogonal, the field is in an incoherent superposition. The 
collapse of Rabi oscillations is a consequence of the de- 
creasing overlap |(a m |a m /)| from Eq. (23). 



Deep in the collapse phase the overlap between differ- 
ent contributions |cr m )(g)|a m ) in Eq. ([8| is negligible. As a 
consequence the phase of ip m (t) in Eq. (25), and thus the 
term At responsible for atomic oscillations, drops out of 
the expressions. The atomic expectation values are now 



(Jz(t))= £ m\(m\0 + 50)\ 2 cos £ m (t) 



(28) 
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FIG. 4. (Color online) Formation of field cat states in the 
CR dynamics for ft/ A = 1/20 <1, with A/A = 0.02, j = 10 
and 0o = 7r/4, cto = 3. Upper panels: Spin expectation val- 
ues (J x / Z ) (left), entanglement entropy S and field variance 
Af (right). Lower panels: Atomic and field Husimi func- 
tion Q a {0,(j)) (left), Qf(a) (middle) and Schmidt coefficients 
(right), at t = 25 x (2n/ft). For Q a (0, 0), we show the hemi- 
sphere < < 7r/2, with = in the center, = ty/2 on 
the outer circle, and < (j) < 2tt running counterclockwise as 
depicted. The gray dashed circles in this and the following 
figures indicate = 7r/4. 



(Jx(t))= £ m\(m\0 + 50}\ 2 sin £m(t) , (29) 



with £ m (t) = SO cos(Q — 2ujEm)t. Expanding cos£ m (£) = 
1 + 0(A 2 ) we see that, up to small corrections, (J z ) in 
the collapse phase is given by the constant value 



(^W)|coiia P se= rn\(m\0 + 50)\< 



m=-j 



= (0 O + S0\J z \0 o + 50) = jcos(# + SO) . (30) 

This result is again a consequence of the rotation of the 
initial atomic coherent state in Eq. ( [l9| ) . For the parame- 
ters in Fig. [l] the expected value is j cos(#o + SO) ~ —0.15, 
which is in full agreement with the numerical results. On 
the other hand, (J x (t)} retains through the sin£ m (£) term 
of order A the cos(^t)-dependence seen in the middle 
panel in Figs. [Tl [3] 



For later times, inspection of Eqs. (22), (23) shows that 
periodic revivals occur at multiples of Te- For Fig. [TJ 
the estimate Te = 100 (27r/fi) is close to the numerical 
value. The appearance of a periodic CR pattern is a 

(2) 

consequence of the linear dependence of Emn on both 
m, n, which is a significant difference from the RWA at 
resonance, where the energy correction ex y'n prevents 
truly periodic revivals [6]. 



IV. GENERATION OF ATOMIC AND FIELD 
CAT STATES 

The non-resonant perturbation theory from the previ- 
ous section shows that in the opposite cases 0«A and 
ft ^> A either the term ex a^aJ z or cx J 2 in the effective 
Hamiltonian H determines the dynamical properties over 
the respectively shorter time scale Te or T5. Both terms 
lead to a different structure of the wave function that 
is related to the generation of either field or atomic cat 
states. These cat states appear as linear combinations of 
well-separated atomic or field coherent states. A "per- 
fect" field cat state would, e.g., be the state \a) ± |— a) 
for \a\ > 1. 

We have discussed the case O < A already in some de- 
tail, and will now revisit the structure of the wave func- 
tion for the situation shown in Fig. [4j As expected, the 
collapse phase coincides with large atom-field entangle- 
ment and large field variance. The coherent state pa- 





FIG. 5. (Color online) Atomic/field Husimi function of the 
two largest contributions in the Schmidt decomposition from 
Fig. [3] with respective weight 0.37 (left), 0.25 (right). 
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FIG. 6. (Color online) Atomic/field Husimi function of the 
two relevant contributions in the Schmidt decomposition at 
t — 60 x (27r/Q) ~ Te/2 (remaining parameters as in Fig.[4|, 
with respective weight 0.65 (left), 0.32 (right). 



FIG. 8. (Color online) Atomic/field Husimi function of the 
two largest contributions in the Schmidt decomposition from 
Fig. [7| with respective weight 0.78 (left), 0.18 (right). 



rameters a m (t) in the wave function Eq. ([8| differ by an 
angle that is a multiple of 27r Whenever t/T E is 
a rational number, some of the a m are equal such that 
fewer coherent field states appear in the wave function. 
This explains the dips in the entanglement entropy S(t) 
in Figs. [TJ [4} 

At t = 25 x (2tt/0) ^T e /5, the field Husimi function 
Qf is a superposition of five Gaussian peaks. For the 
atomic state, a similar fivefold pattern cannot be clearly 
identified in the atomic Husimi function [15] 



(31) 



Q a (e,<j>) = \(e,<j>\i>)f 



which is defined via atomic coherent states from Eq. ( |17| ) 
similar to the field Husimi function. It gives the pseudo- 
spin probability distribution as a function of spherical 
angle coordinates 0, <p. For A = 0, the phase space sphere 
would rotate rigidly around the z-axis passing through 
the origin = 0. 

If the atomic states \a m (t)) in Eq. ([8| were mutually 
orthogonal the field superposition seen in Qf would be 
completely incoherent. We know from perturbation the- 
ory that the states are not orthogonal due to the rotation 
with E/t in Eq. (24), which allows for a coherent superpo- 



sition and the appearance of field cat states. To check up 
on this possibility we use the Schmidt decomposition of 
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FIG. 7. (Color online) Formation of atomic cat states for 
n/A = 20 > 1, with A/A 0.5, j 10 and = tt/4, a = 3 
in the initial state. The panels show the same quantities as 
in Fig. [3] with t = 32 x (2n/A) in the lower row. 



the atomic/field wave function [17]. The Schmidt coeffi- 
cients, i.e. the eigenvalues of both the atomic/field den- 
sity matrix depicted by Q a /f-> have five relevant contri- 
butions. The two largest are shown separately in Fig. [5j 
We can identify a field cat state in the largest contribu- 
tion (left panels), while the corresponding atomic state 
is strongly squeezed but not a cat state. In the second 
largest contribution the indication of an atomic cat state 
is visible. We also conclude that the fivefold field super- 
position in Fig. [4] is partially coherent. 

Halfway through the collapse phase, for t « T#/2, 
the two contributions shown in Fig. [6] comprise 98% of 
the wave function. Now, atomic and field cat states 
appear simultaneously. Note that the appearance of 
field cat states is again a consequence of the transfor- 
mation W in the perturbative result and the resulting 
non-orthogonality of the \a m ) states in Eq. Other- 
wise, the wave function would have the form~[cr+)|a) + 
|cr_)|— a) with two orthogonal atomic states |cr±), and no 
field cat states could appear. Instead, the final rotation 
with W in Eq. (24) leads to a finite overlap (cr + |a_) ^ 0, 
and field cat states ~ \a) ± |— a) occur in the Schmidt 
decomposition. 

In the opposite case ft ^> A, addressed in Fig. [7| the 
wave function remains close to a product form \<r(t)) (8) 
\ot(t)) for times t < Ts <C T#. The entanglement entropy 
and field variance are correspondingly small. In contrast 
to the initial short-time dynamics of Rabi oscillations in 
Eq. ( 26 ) that we deduced from the product form for times 
t <C Ts,Te, the term J% in H now gives rise to squeezing 
and splitting of the initial atomic coherent state. In the 
situation shown in Fig. [7| the atomic Q-function at t = 
32 x (2tt/A) w T s /2 <C T e has a two-fold structure. It 
appears since, in the present case of integer even j, 



exp 



iU 2 

2 z 



\0) 



l\0) 



(32) 



is a linear combination of two atomic coherent states 
\—0). The two relevant contributions in the Schmidt de- 
composition are shown in Fig. [8] As opposed to the previ- 
ous case Q <C A, we identify the signatures of an atomic 
cat state together with a single coherent field state in 
the largest contribution (left panels). The annulus in Qf 
for the second largest contribution (right panels) origi- 
nates from a superposition of two field coherent states 
\a ± 5a) with \Sa\ <C 1. This structure is a precursor of 
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0.5 

t / (2%/Q) 



12 3 4 

t / (2iUQ) 



FIG. 9. (Color online) Comparison of the SCA to the exact 
dynamics of the Rabi model (j = 1/2) with Q/A = 0.01, 
A/A = 0.02 and 0o = tt/2, ao = 5. The left panels show 
(J z (t)) over a short time interval, presenting the exact nu- 
merical, perturbative and SCA result together with the result 
in the limit of a field-driven atom (Eq. (33)). The right pan- 
els cover a longer time interval and include the entanglement 
entropy S(t) (red curve) from the numerical calculation. 

the field state splitting through the a^aJ z -teim for later 
times T s <t< T E . 



V. ATOM-FIELD ENTANGLEMENT AND THE 
SEMI-CLASSICAL APPROXIMATION 

The CR mechanism for Q <C A does not depend on the 
non-classical properties of a quantized field, but only on 
the possibility of atom- field entanglement. We thus ex- 
pect that the basic signature of this mechanism, the col- 
lapse of Rabi oscillations, occurs also for atoms in strong 
radiation fields close to the classical field limit where the 
field quantization plays no role. 

The classical field limit can be defined rigorously as the 
limit |ao| — > oo, keeping X\ao\ constant. In this limit, 
the field mode evolves independently of the atoms be- 
cause the product X(J X ), which gives the strength of the 
atomic influence on the field, goes to zero for A — >• and 
finite j. The product A (a + a^), which determines the 
influence of the field on the atoms, remains finite. The 
Dicke model reduces to the model of an atomic ensemble 
driven by an external field B(t) = (2Aao cos Qt, 0, — A). 
The atomic expectation values J = (( J x ), ( J y } : (J z )) obey 
the equation of motion 



d t J = B(t) x J . 



(33) 



The characteristic signatures of such a driven system are 
Rabi oscillations. 

For finite a, A the field mode does not evolve indepen- 
dently of the atoms, but we expect only small corrections 
from the classical field dynamics for small A. One at- 
tempt to include these corrections is the semi-classical ap- 
proximation (SCA) (see, e.g., Ref. [20 for a discussion). 
The SCA is based on the assumption that the coupled 
atom- field system remains in a product (j)(t))®\a(t)) 
of coherent states during time evolution. This assump- 
tion allows for decoupling of the equations of motion for 
the atomic and field expectation values. 

After decoupling, the atomic state evolves again simi- 
lar to a spin in a magnetic field B(a) = (2ARea,0, — A). 
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FIG. 10. (Color online) Comparison of the SCA to the exact 
dynamics, still with Q/A = 0.01 and 0o — tt/2 as in the 
previous Fig. [9] Left panels: Results for the Rabi model (j = 
1/2) with A/ A = 0.005 and a = 20. Right panels: Results 
for the Dicke model with j = 10, and other parameters A/A = 
0.02 and O = tt/2, a = 5 identical to Fig. [9] 



Now, however, the field state evolves as for an oscillator 
with an additional external force X(J X ) that accounts for 
the back-reaction from the atomic ensemble. The corre- 
sponding SCA equations of motion are 

9 t J = B(ft)xJ, id t a = {la + \(J x ) . (34) 

The SCA equations of motion become exact in the clas- 
sical field limit, where they reduce to Eq. (33). 



From the equations of motions we can observe a poten- 
tial problem of the SCA that arises from the generation 
of atom- field entanglement in the true Dicke dynamics. 
The change of the field state in Eq. ( |34| ) is oc A, and the 
influence on the atoms ex A 2 . Therefore, the dynamically 
relevant time scale in SCA is cx 1/A 2 , just as the scaling 
of the entangling time Te from Eq. (14). We should thus 



expect that for time scales on which the SCA differs from 
the simpler Eq. (33) significant atom- field entanglement 



has been generated for which the SCA cannot account. 

We consider exemplarily the limit Q — >• of a classi- 
cal field with negligible energy quantization. In Fig. [9] 
we compare the SCA with the exact and perturbative re- 
sult for the dynamics, and with the simplified Eq. ( [33] ) of 
a driven atomic ensemble. For short times (left panels) 
all four descriptions agree and describe Rabi oscillations, 
which are characteristic for classical field dynamics. For 
longer times (right panels), we see that significant entan- 
glement is generated over the first few field periods even 
in the most simple Rabi case j = 1/2. The SCA cannot 
account for entanglement and accordingly misses the col- 
lapse of Rabi oscillations entirely. We note that the SCA 
simply reproduces the result also obtained with the sim- 
pler Eq. ( 33 ) . Note also that we are in the weak coupling 



regime, where the non-resonant perturbation theory de- 
scribes the dynamics accurately. The failure of the SCA 
is not a result of strong atom-field coupling. 

A situation with larger ao is shown in the left panel 
in Fig. [To] The collapse of Rabi oscillations takes place 
for later times since Te is larger, but both SCA and the 
simpler Eq. ([33] ) again fail in the same way. In the right 
panel in Fig. [Tojwe see that this behavior is not restricted 
to the Rabi case, but occurs equally for j ^> 1/2. We note 
that its significance is not reduced since the present vio- 
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lation of classical field dynamics is not a consequence of 
quantum fluctuations of the atomic system, which would 
decay with The relevance of entanglement, bounded 
by ln(2j + 1), even increases with j. 

The violation of the mean-field assumption of negli- 
gible atom-field correlations already at weak coupling 
presents a genuine problem for the SCA. As a conse- 
quence of entanglement generation, the SCA can be cor- 
rect only as long as the field dynamics remains decoupled 
from the dynamics of the atomic ensemble. In that sit- 
uation the atomic ensemble is already described by the 
simpler Eq. (33): The SCA does not improve on a model 



with entirely classical field dynamics where any influence 
of the atoms on the field state is neglected. We conclude 
that the SCA fails to account for the back-reaction of the 
atoms on the field, even at weak coupling or large j. 

The non-resonant CR dynamics discussed here is an 
example of non-trivial coupled quantum-classical dynam- 
ics [21] , as is evident from the wave function Eq. ([8| with 
multiple coherent field states. In our case, the signatures 
of coupled quantum-classical dynamics are the periodic 
CR patterns which cannot be explained in a simple mean- 
field description. 



VI. CONCLUSIONS 

In summary, our analysis explains the CR patterns 
of Rabi oscillations in the non-resonant Dicke model by 
the dynamical splitting of atomic or field coherent states, 
which takes place on different time scales distinguished 
by the detuning ratio Q/A. The collapse phase is ac- 
companied by the buildup of atom-field entanglement or 
atomic squeezing, whose subsequent decay leads to peri- 
odic revivals. Both atomic and field cat states can arise 
in the collapse phase. The quantitative explanation of 
this behavior is provided by perturbation theory for large 
detuning. Application of the RWA to the non-resonant 
case would erroneously predict a single time scale (cf. 
App. |A| ), instead of the two time scales obtained in the 
correct calculation. 

The non-resonant CR patterns arise through a dynam- 
ical mechanism that involves the creation of highly non- 
classical states from initial classical state preparations. 
They give direct evidence for quantum entanglement and 
coherent quantum superpositions in atom-field and re- 
lated systems. Even close to the classical field limit, 
where other quantum properties such as field quantiza- 
tion are of minor importance, atom-field entanglement 
can prevail over the semi-classical dynamics that would 
occur for a hypothetical non-correlated atom-field sys- 
tem. This indicates how semi-classical approximations 
can fail because they neglect atom-field correlations. 

The observation of non-resonant CR patterns might 
become possible in experiments using superconducting 
circuits instead of optical cavities [22]. These experi- 
ments can reach the regime of strong coupling or large 
detuning [23 , which requires calculations beyond the 



RWA [24] [25] . Although the experimentally controllable 
detuning can be made large, the most serious obstruction 
against observation of the dynamical patterns described 
here is the necessity of preserving quantum coherence 
over many Rabi oscillations. Further improvement of ex- 
perimental techniques may resolve this issue. 
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Appendix A: Non-resonant perturbation theory for 
Dicke-type Hamiltonians 

We give here the result of second order perturbation 
theory for Dicke-type Hamiltonians 



H = ujJ z + Qafa 

+ g(aJ+ + a) J_) + g(a) J+ + aJ_ 



(Al) 



with J± = J x ± iJ y denoting spin ladder operators, in 
the non-resonant case \uj\ ^ Q. 

Standard perturbation theory gives a correction to the 
eigenstates and eigenvalues of the non-interacting Hamil- 
tonian in the first line in Eq. (Al). The result can be 
expressed in the form of Eq. (10), with a unitary trans- 
formation 

Z7 = exp[T!+T 2 ] (A2) 
that accounts for the change of the eigenstates, where 
9 , J + 7 n ■ 9 



Ti = 



T 2 



-{aJ + - at J_) + -^(atj+ - aJ_) , (A3) 



99 



n 2 



2uj 



(Jl - Jl 



(at 2 



(A4) 



and an effective Hamiltonian 

H = wJ z + QaJa 

g 2 (uj + fi) + g 2 (co - n) 



n 2 



g 2 (u + n)-g 2 (uj-n) 

oj 2 -n 2 



(2a f a + l)J z 
{J 2 ~ Jl) 



(A5) 



that gives the perturbed eigenvalues. By construction, 
H is diagonal in the unperturbed eigenstates and con- 
tains only operators J z and a^a, and the constant of mo- 
tion J 2 . Perturbation theory thus provides us with an 
approximate diagonalization of the Hamiltonian, in the 
sense that both sides of 



U^HU = H + 0({g,g} 3 ) 



(A6) 



differ by terms of third or higher order in the coupling 
constants. Alternatively, the contribution of T 2 can be 
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included in the effective Hamiltonian as in Ref. [25 J. The 
present formulation with a diagonal H is preferential for 
the study of the dynamical evolution of the wave function 
since it allows for direct evaluation of Eq. (10). 



Our Hamiltonian from Eq. ([I]) corresponds to the 
choice g = g, including counter-rotating terms in H. 
The second order term T2 is finite and contributes to the 
atomic and field squeezing, but it appears as a higher or- 
der correction to the leading order term T\ in U . There- 
fore, it is not relevant for the understanding of the dy- 
namical effect discussed in the main part of the text and 
was not included in Eq. (12). Within the perturbative 



setting the relevant second order term in H, which grows 
during time evolution, is separated from the irrelevant 



second order term in U . 

In the RWA, applicable at resonance lj « O, counter- 
rotating terms are dropped from the Hamiltonian by set- 
ting g = 0. In Refs. [I3j [14] an even simpler model with 
g = and ft = was taken as the starting point. Within 
RWA, the term T2 is zero in accordance with the fact that 
the number of excitations (corresponding to the operator 
J z + a^a) is conserved. More importantly, the effective 
Hamiltonian now contains the same prefactor g 2 / (uj — Q) 
in front of the two operators a^aJ z and J|. Instead of the 
two different time scales for atomic and field squeezing 
introduced in Eq. ( 14 ) only a single time scale appears in 



the RWA. By construction, the RWA is valid only close to 
resonance and incompatible with the non-resonant per- 
turbation theory. 
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